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Abstract 

Equilibrium particle densities near a hard wall are studied for 
a quantum fluid made of point charges which interact via Coulomb 
potential without any regularization. In the framework of the grand- 
canonical ensemble, we use an equivalence with a classical system of 
loops with random shapes, based on the Feynman-Kac path-integral 
representation of the quantum Gibbs factor. After systematic resum- 
mations of Coulomb divergences in the Mayer fugacity expansions of 
loop densities, there appears a screened potential <f>. It obeys an 
inhomogeneous Debye-Hiickel equation with an effective screening 
length which depends on the distance from the wall. The formal 
solution for <J> can be expanded in powers of the ratios of the de 
Broglie thermal wavelengths X a 's of each species a and the limit 
of the screening length far away from the wall. In a regime of low 
degeneracy and weak coupling, exact analytical density profiles are 
calculated at first order in two independent parameters. Because 
of the vanishing of wave-functions close to the wall, density profiles 
vanish gaussianly fast in the vicinity of the wall over distances Ac's, 
with an essential singularity in Planck constant h. When species 
have different masses, this effect is equivalent to the appearance of a 
quantum surface charge localized on the wall and proportional to h 
at leading order. Then, density profiles, as well as the electrostatic 
potential drop created by the charge-density profile, also involve a 
term linear in h and which decays exponentially fast over the clas- 
sical Debye screening length £ D . The corresponding contribution to 
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the global surface charge exactly compensates the charge in the very 
vicinity of the surface, so that the net electric field vanishes in the 
bulk, as it should. 

KEYWORDS : Coulomb interactions, quantum mechanics, 
hard wall, grand-canonical ensemble, inhomogeneous Debye equa- 
tion, surface charge. 



1 Introduction 

1.1 Issue at stake 

In the present paper the equilibrium density profiles in a quantum fluid of 
point charges are studied in the vicinity of an impenetrable hard wall. The 
wall, which occupies the semi-infinite region x < 0, has no internal structure 
and its dielectric constant is the same as that of the medium where charges 
move. On the contrary, the fluid made of n a particle species is described at 
the microscopic level in the framework of quantum statistical mechanics. 
Two point charges e a and e a > (where a is a species index) interact via the 
electrostatic interaction e a e a > v(r — r'), where 

v ^-^ = ]^\ (L1) 

in Gauss units. (The charge e a includes a factor l/^fe^ in energy terms 
when charges are embedded in a continuous medium with a relative dielec- 
tric constant e m with respect to the vacuum.) The interaction is transla- 
tionally invariant, and the anisotropy lies only in the geometric constraint 
enforced by the presence of the wall. The exact analytical expressions of 
the density profiles p a (x) 's are obtained in a regime where exchange effects 
are negligible and where Coulomb coupling is weak. Results hold for the 
electron-hole gas in an intrinsic semi-conductor in the vicinity of a junction 
or for a dilute and hot quantum plasma near a vessel wall. 

The interesting point of the model is that it exhibits how a quantum 
charge effect, gaussianly localized over de Broglie thermal wavelengths in 
the vicinity of the wall, is carried by long-range Coulomb interactions up to 
larger distances from the wall, with an exponential decay over a scale equal 
to the coulombic screening length. Indeed, whereas the density in a classical 
ideal gas is uniform in the whole region x > and is discontinuous on the 
wall surface, the quantum density is continuous and vanishes at x = 0, 
because of the continuity of wave-functions and their cancellation inside 
the impenetrable wall. At the inverse temperature f3 = l/k B T (where k B 
is Boltzmann constant), in a low-degeneracy limit quantum statistics is 
reduced to Maxwell-Boltzmann statistics, and the density (x) of species 
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a in an ideal gas with quantum dynamics vanishes gaussianly fast over the 
scale of the thermal de Broglie wavelength A a , 

=P°(l- e^ 2 ^) + plO (j^fj ■ (1-2) 
In II1.2II is the bulk density for species a, and 

A Q = tuf^-, (1.3) 



where h is Planck constant and m a is the mass of species a. The quantum 
expression with Maxwell-Boltzmann statistics 11.211 is valid up to terms 
of order (A a /a Q ) 3 , where a a is the mean interparticle distance between 
particles of the same species a ((4/3)7rp^a 3 = 1). 

When Coulomb interactions are taken into account, there arises only one 
third typical length scale, because Coulomb interactions are scale-invariant. 
Then only two independent dimensionless parameters rule the physical 
regimes of the system: the degeneracy parameter and the Coulomb cou- 
pling parameter. The degeneracy parameter is (A/a) 3 , where A = sup a {A Q } 
and a is a typical mean interparticle distance. When exchange effects are 
negligible, 

J «1. (1.4) 

The third length scale may be chosen to be either the classical screening 
length £ D 



\ 
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4^]Tpge 2 (1-5) 

a=l 



or the Landau length /3e 2 , namely the classical closest approach distance 
between two typical like-charges e with kinetic energy of order 1/ 13. Hence- 
forth, the classical Coulomb coupling parameter can be chosen to be equal 
either to the ratios (a/^ D ) 3 , /3e 2 /a = T, or (3e 2 /£_ D = 2e D . These ratios are 
proportional to one another, 

f]\e D o,Y^. (1.6) 

More precisely they are linked by the relations (a/£ D ) 3 = Ce D and e D = 
[C/8] 1 / 2 T 3 / 2 , where C is a numerical factor which depends on the com- 
position of the fluid through £d and on the relation between a and bulk 
densities. In a weak-coupling regime 

£)*«!. (1.7, 
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When both 11.41 and l|1.7fl are satisfied, 

A<a<£ D . (1.8) 

Therefore, the low-degeneracy and weak-coupling regime is also a regime 
where k d X <C 1. 

1.2 Results 

In the low-degeneracy and weak-coupling regime defined by ljl.41) and l|1.7l) , 
the analytical expression for the profile density p a (x) is calculated in a 
subregime where the first coupling correction, of order e D — K D /3e 2 /2, and 
the first diffraction correction, of order k d A, dominate other coupling and 
exchange corrections. At order e D classical contributions do not involve the 
short-range cut-off that must be introduced in order to prevent the collapse 
of the system in the limit where h tends to zero. As shown in Section El 
the subregime corresponds to a scaling where 

el < (-) «e d . (1.9) 



a J 

tQ} can be reexpressed as s 3 D < (k d X) 3 <C e 2 D , since (A/a) oc K D A/e]/ 3 
by virtue of 11. fill . In the subregime tP} the density profile in the region 



x > reads 



+ p B a o{s d ,k d \) (1.10) 



where o(e Dl n D \) denotes a sum of terms which tend to zero faster than 
either e D or k d X when these parameters vanish (See 115. 21 ). When the 
latter terms are neglected, p a (x) appears as the product of the ideal-gas 
density 111.211 with a function arising from interaction corrections. (The 
generic properties of density profiles are discussed in Section H3) 

The density profile l|1.10|) results from the combination of three effects: 
first, the vanishing of quantum wave-functions in the vicinity of the wall; 
second, the geometric repulsion from the wall QJ, described by the classical 
part of the screened self-energy due to the deformation of screening clouds, 
(l/2)n D /3e 2 t L(u), given in US. 141) : third, the interaction e Q $(x) with the 
electrostatic potential drop $(a:) with respect to the bulk, which is created 
by the charge density profile J2-y e iP-y( x ) itself. (The sign of the latter 
interaction depends on the sign of e a .) The potential drop $(x) in 11.101 
is the sum of a classical contribution and a quantum "diffraction" effect, 
linear in H, 

$(x) = $ c1(£d) (x) + $q u ^ A ) ( x ). (1.11) 
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$ci(e. D )( a ;) ) Q £ orc i er eD /(^ e ) ; i s written in H5.40I) and ^^'^^(j;), of order 
k d X/ (/3e), reads 

$qu( KD A) (a;) = _^ e -^ withg ^ _^S 7 (e 7 /^)P? ( } 

( j,qu(K D A)j- 2 ,^ appears only when species have different masses, because of 
the bulk local neutrality 

5>aP2=0. (1.13) 

a 

The density profile ll.lfljl can be rewritten as 

Pa (x) = (l-e- 2x2 l^) [p*")(i) + We a Be- toI ] + p B a o(e D , Kn X), 

(L14) 

where /0a (a;) is the classical density profile calculated up to relative or- 
der e D pQ and written in Q5.47I1 . We stress that the direct contribution 
ljl.2ll from the vanishing of wave-functions in the ranges A Q 's from the wall 
has an essential singularity in h, whereas the quantum part of the elec- 
trostatic potential at leading order, $q u ( K o A )(a;), is linear in H. (We recall 
that for systems invariant under translations - which is not the case here - 
and with sufficiently smooth potentials - such as the Coulomb interaction 
- ^-expansions involve only even powers of h, as can be seen for instance 
in Wigner-Kirkwood expansions.) 

The appearance in the electrostatic potential of a h-term which de- 
cays exponentially fast over the classical Debye screening length £ D has the 
following physical interpretation. When species have different masses, the 
global charge tr< carried by the fluid (per unit area) over all distances x < a 
from the wall is essentially created at leading order by the differences in the 
Gaussian density profiles and is concentrated over a width of order A. Since 
A is negligible with respect to the bulk mean interparticle distance a, the 
leading-order charge a ^( KDX > can b e seen as a surface charge localized at 
x = 0. As shown in Section IQ1 the surface charge a ^ KDX ^ which appears 
even in the zero-coupling limit, creates an electrostatic potential through 
the classically-screened Coulomb interaction (calculated at leading order), 
and this potential is equal to the leading ft-term $i u ( KdA ) (x) in the elec- 
trostatic potential created by the charge-density profile ^ Q e a p a (x). 
Moreover $i u (' tDA ) (a;) is involved in the density profiles in such a way that 
the leading-order global charge (7 ^ u ( KdA ) carried by the fluid (per unit area) 
over all distances x > a, and which is dilute over the scale £ D , compensates 
^qufKcA) ^ gee Section lfi^4l . Indeed, since the wall is made of an insulating 
material and carries no external charge, the global surface charge a carried 
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by the fluid per unit area vanishes at equilibrium |2J 



oo 




(1.15) 



^qufnuA) an( j (jjquffiDA)^ _ g-j are es ti m ated in the case of the intrinsic 
semiconductor GaSb. The case where the wall has not the same dielectric 
constant as the medium where the fluid is embedded is commented in 
Section □ 

1.3 Methods 

Before going into details, we summarize the general methods displayed in 
Sections EH3 

First, a formalism based on path integrals and devised for the study of 
bulk properties in Coulomb fluids - with Maxwell-Boltzmann statistics 
then quantum statistics is generalized to a semi-infinite geometry (Sec- 
tion EJ, The system is studied in the grand-canonical ensemble (Section 
12. li . A degeneracy of physical quantities with respect to fugacities arises 
from the neutrality constraints enforced by the long-range of Coulomb in- 
teractions. We investigate the nature of this degeneracy, and we show that 
we are allowed to split the latter degeneracy in order to impose the local 
neutrality in the zero-coupling limit fSection I2.2JI . (This trick allows to 
simplify weak-coupling expansions performed in Sectional) By use of the 
Feynman-Kac formula ( Section 12.311 . quantum dynamics can be described 
by a functional integral over Brownian paths, which correspond to quan- 
tum position fluctuations. As in the bulk situation, the quantum system 
of point charges is equivalent to a classical system of loops with random 
shapes (Section 12. 4|) , The only difference in formulae for the bulk or for 
the vicinity of the wall is that the path measure is anisotropic and depends 
on the distance from the wall in the second case. 

Then methods originally devised for classical fluids with internal degrees 
of freedom can be used (Section EJ. In Section l3~Tl we introduce general- 
ized Mayer diagrams for the fugacity expansion of the loop density of each 
species. Point weights in those diagrams depend both on the internal de- 
grees of freedom of loops - charge and shape - and on the distance x from 
the wall. Because of the long range of Coulomb interaction, every Mayer 
diagram that is not sufficiently connected corresponds to a divergent inte- 
gral in the thermodynamical limit. These divergences disappear after exact 
systematic resummations analogous to those performed in Ref.@] (Section 
13.211 . (Details are provided in Appendix El) Resummations introduce a 
screened potential <j>(r, r'), solution of an inhomogeneous Debye equation 




k 2 (x)] 0(r, r') = -47T<5(r - r'), 



(1.16) 
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where the effective screening length 1/k(x) depends on the distance x to 
the wall because of the vanishing of wave-functions at the wall surface (see 

CE2J)- 

At this point the difficulty to be circumvented is the resolution of equa- 
tion Ijl.lfiH (Section I3J. The equation can be turned into a one-dimensional 
differential equation by considering the Fourier transform 
4>{x, x' , k) of <p{x, x' , y) in the directions parallel to the wall surface (Section 
14. 11) . Let (j)^ ( x , x \ y) be the expression that $(x,x',y) would take if the 
profile k(x) were uniform and equal to its bulk value k in the region x > 0. 
A formal series representation of the solution 4>(x, x' , k) — tfi^ (x, x' , k) has 
been given in Ref. pj, where a similar equation arises in the case of a classi- 
cal charge fluid in the vicinity of a wall with an electrostatic response. This 
series provides an expansion of the solution (f>(x, x' , k) around <^°) (x, x' , k) 
in powers of the small parameter kX, where k is the limit of k(x) when x 
goes to infinity, while A is the length scale over which k(x) varies quickly 
when x approaches 0. The expansion of <p(x, x' : k) — </>(°) (x, a;', k) in powers 
of kA is uniform in x and x' . 

In the low-degeneracy and weak-coupling regime to be studied (Section 
15 .111 , the condition ft A <C 1 is met. The screened self-energy is purely 
classical at leading order and the quantum correction appears only at order 
ex ftA, where e is defined as e D with k in place of k d , e = (l/2)K/3e 2 (Section 
15.21) . A scaling analysis performed in the low-degeneracy and weak-coupling 
limit (Section and Appendix [Bj shows that only one resummed Mayer 
diagram contributes to density profiles at first order in e and kX. The 
electrostatic potential drop &(x) created by the charge-density profile is 
identified in the formal expression of the contribution from this diagram 
fSection I5.4JI . Because of the local neutrality condition in the bulk, only 
the classical zeroth-order term in the ftA-expansion of 4> — 4>^ proves to 
contribute to the potential drop $(x) at leading orders e and kX. 

2 General formalism 

2.1 Grand-canonical ensemble and statistics 

We recall that we consider a fluid made of n s species (indexed by a) , each 
of which is characterized by its mass m a , its charge e a , and its spin S a h. 
(In the following, interactions involving spins will be neglected and spin 
will only determine the nature of quantum statistics.) The Hamiltonian 
operator H{N a } of a system which contains N a particles of each species a 
reads 
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({N a } is a shorthand notation for {N a } a= i j ... t n s and the particle index i 
runs from 1 to N — J2 a ^ a -) ^ n Q2.1fl the first term which involves the 
momentum operator p is the global kinetic energy of the system. The 
second term is a sum of one-body potentials V SR (xi) which describe the 
particle- wall interactions. We choose a simple classical hard-wall modeliza- 
tion, where the atomic structure of the wall is ignored. The effect of the 
wall is only to prevent particle wave-functions from propagating inside the 
negative-x region occupied by the wall, 

The wall repulsion is independent of the particle species. The sum of pair 
interactions in the third term involves only Coulomb potential 1)1. 11 . 

The fixed parameters of the system are the volume |A| of the region 
A that the fluid occupies, the area \d w A\ of the fluid-wall interface, the 
temperature, and the densities p%'s far away from the boundaries of A. 
We use the grand-canonical ensemble where the parameters are the volume 
|A|, the area |<9 W A|, the inverse temperature (3, and the chemical poten- 
tials {fJ- a (V R )} a =i,...,n s of particles in a reservoir where the electrostatic 
potential takes the uniform value V R . The grand partition function reads 

SOS, W,|A|,|a w A|) = £ T^je-'fr^-E-^H (2.3) 

{Nc} 

where the trace Tr^™^ , is restricted to the quantum states that are prop- 
erly symmetrized according_to the Bose-Einstein or Fermi-Dirac statistics 
obeyed by each species. (N a is the particle-number operator for species 
a.) As in the classical case, the density profile p a (x) can be determined 
by using a functional derivation of S [ju^], where ^ Q p a N a is replaced by 
f A drji^(x)p~^(x). The relation is 



o (x)- hm 1 SlnE ^ 

U \ / 11111 



|A|->+oo (3 SfJ, a (x) 



(2.4) 



The formalism and results presented in SectionsElandElcan be obtained 
with quantum statistics (as detailed in Section l2~4*jl . However, our aim is 
to produce explicit analytical results in the low-degeneracy regime l|1.4ll . 
We have checked that quantum statistics effects arise only at order (A/a) 3 . 
In other words [Hj 

3 = Hmb+O^Q) , (2.5) 

where the Maxwell-Boltzmann grand partition function Smb is a trace over 
tensorial products of one-particle wavefunctions which are not symmetrized 
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according to species statistics. If the tensorial product of N = J2 a N a one- 
particle states in position representation is denoted by |{r.j}), the grand 
partition function Smb, where only dynamics is quantum, reads 



SmbGs.W.Iau^ad 



n 



\JV„ 



NJ. 



N 



Y[drt 



(2.6) 

({rj|e-^^>|{rj). 



where 2S a + 1 is the spin degeneracy factor, which arises because spin does 
not appear in the expression 112. If) of the Hamiltonian. In l|2.6l) we have 
used the commutativity of the operators -H{iv a } and N a 's. 



2.2 Degeneracy with respect to fugacities 

In the following we will take advantage of a degeneracy of physical quanti- 
ties with respect to fugacities that arises from the vanishing of the global 
volumic and surfacic charges of the system in the thermodynamic limit. If 
(• • • ) denotes a grand-canonical average, the thermodynamic limit of the 
charge in the fluid is 

lim(^ e a N a ) = ( e *Pa ) I A| + a\d w A\ + o (\d w A\) , (2.7) 

where o(|<9 w A|) denotes a term which diverges more slowly than the area 
|9 W A| when the size of the domain A becomes infinite. The expression 
of a in terms of the thermodynamic limits of density profiles is given in 
Q1.15J1 . As a consequence of the existence of the thermodynamical limit [gj, 
the macroscopic volumic charge (J2 a e aPa) \M vanishes, and, in the case of 
an insulating hard wall that is not externally charged, the surfacic charge 
in the fluid <r|<9 w A| is also equal to zero (whether the dielectric constants 
in the wall and in the medium where the fluid is embedded are equal or 
not). Indeed, in the grand canonical ensemble |12..'>L the summation over 
microscopic states involves non-neutral configurations, but the self-energies 
of these globally charged configurations give them exponentially vanishing 
weights in the thermodynamic limit, because they are not compensated by 
interaction energies with external charges inside the walls. 

Since the bulk charge neutrality is satisfied for any set of chemical 

potentials, the bulk densities p^'s are determined by only n s — 1 indepen- 
dent functions of the n s chemical potentials ^ Q 's. In the present paragraph 
we investigate more precisely the nature of the corresponding degeneracy. 

By definition, the electrostatic energy in the Hamiltonian H2.1H used in 
5 H2.3H is the difference between the electrostatic energy of the interacting 
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system and the energy V R ^2 a e a N a of the noninteracting system in the 
reservoir where the electrostatic potential takes the uniform value V R . In 
other words, the dependence of chemical potentials with respect to the 
potential V R in the reservoir is just 

» a (V R ) = fi a (0) + e a V R . (2.8) 

The global volumic and surfacic neutralities are linked to the invariance 
of the thermodynamic limits of observables under a translation of the origin 
for the electrostatic potentials. Indeed, if the latter origin is translated by 
an amount —AV, then, the reference potential V R of the reservoir becomes 
V R + AV, the Hamiltonian is unchanged (since the insulating wall carries 
no external charge) , and the only change in 3 H2.3J1 is an extra contribution 
AVJ2 a e aN a arising from J2 a Va(V R )N a . Then the thermodynamic limit 
of In 3 is increased by 



Urn In 3 1 = (3AV hm(^ e a N a ) (2.9) 



According to H2.7I) . (I1.13|) . and II1.15J1 . the latter variation vanishes up to 
order \d w A\ included. 

Since a translation — AV^ of the origin for the electrostatic potentials is 
equivalent to an increase e a AV of every /i Q , the nature of the degeneracy 
of physical quantities with respect to chemical potentials is that physical 
quantities are invariant under the addition of an energy e a AV to every 
chemical potential fi a . The corresponding degeneracy with respect to fu- 
gacities z a 's comes from the definition 

(2S + 1) 

z a (V R ) ee a ' exp W a {V R )\ , (2.10) 
(2ttA£) 

The dependence of fugacities upon V R is given by H2.8J1 . The system in- 
volves charges of both signs, so that the continuous function f(V R ) = 
J2 a e aZ a (V a ) varies from — oo up to +00 when V R varies from — oo to +oo. 
Therefore there exists a value of V R which fulfills the condition f(V R ) = 0. 

As a consequence, since physical quantities are invariant under a trans- 
lation of V R in the fugacities, we can choose a set of fugacities which ensures 
that the local charge neutrality in the bulk is enforced even in the zero- 
coupling limit, namely we can arbitrarily split the degeneracy with respect 
to fugacities by imposing 

^e Q z Q = 0. (2.11) 

We notice that, as shown in Ref.0, in the case of an insulating wall with an 
external charge or in the case of a conducting wall, which becomes charged 
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by influence, the global neutrality of the full system (the fluid plus the 
wall) in the thermodynamic limit implies that condition p.llfl can also 
be fulfilled when S is written with the full Hamiltonian. The "neutrality" 
condition about fugacities <I2 . 1 If) will cause major simplifications in the 
following calculations. 



2.3 Feynman-Kac formula 

The non-commutativity between the kinetic and interaction operators in 
H2.6J1 can be circumvented by using Feynman-Kac formula (HI El- The quan- 
tum Gibbs factor can be rewritten in terms of path integrals, 



x exp 



IM> = 



n 



i 



-(3^2e az e ai / dsv (ri + \ ai £t(s) -rj - A^^s)) 
, Jo 



l<3 



(2.12) 



The kinetic part of the Hamiltonian and the particle-wall interaction are 
taken into account in the measure T> XijClf (£J of the closed Brownian path £ 4 
with dimensionless abscissa s: £(s = 0) = £(,s = 1) = 0. The random path 
A Qi £j(s) with typical extent \ ai describes the quantum position fluctua- 
tions of particle i at position r^. As discussed in Section l2~4l the interaction 
between paths on the r.h.s. of J2.12I is not the usual Coulomb interaction 
between charged wires, since it involves only path elements with the same 
abscissa s. 

The repulsion from the wall causes the anisotropy of the Brownian- 
path measure. The constraint about the quantum particle position, which 
is described by V SR (x), enforces that the x-component £ x of vector £ obeys 
the inequality 

x + \ a £ x (s)>0 (2.13) 
for every s between and 1. The Brownian path measure can be factorized 



(2.14) 



where £m is the projection of £ onto the wall. As in the bulk, the Gaussian 
measure is independent of the position r. Moreover it is rotational- 

invariant and normalized to unity 



/2>(€||) = 1. 



(2.15) 
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On the contrary, as recalled in Ref. ^01; the measure T> Xj0l (£ x ) depends on 
x. with for instance 



= 1 - e-* 8 /A «. (2.16) 
Moreover, the mean extent of the path in the x-direction does not vanish 

J X>«,a(&) & + "• (2-17) 

All moments of the measure tend gaussianly fast to their bulk values over 
the scale of the de Broglie wavelengths. For instance 



ds / X> x ,a(6z) £r(s) 



7T / X 



- v -V, / ' \ A a y 



(2.18) 



where Erfc(u) is the complementary error function defined as 

2 f 00 2 
Erfc(u) = —= / dte~ t2 . (2.19) 
V 71 " A 

Erfc(u) decays as exp[— M 2 ]/(uy / 7r) when u goes to +oo. 

2.4 Equivalence with a classical gas of loops 

In the present paragraph we recall that the quantum grand partition func- 
tion S for point particles can be written as a classical grand partition 
function Si oop for randomly shaped loops @|. The latter formalism in- 
cluding quantum statistics allows one to retrieve property l|2.5|l : quantum 
statistics effects appear only at order (A/a) 3 in the low-degeneracy regime. 
In other words, results at leading order in the low-degeneracy regime lll.4f) 
are the same when the starting partition function is written with Maxwell- 
Boltzmann statistics. Therefore, for the sake of simplicity, we shall directly 
consider Smb and we shall drop the index MB from now on. 

The quantum grand partition function Il2.(ili can be rewritten by use of 
the Feynman-Kac formula 112. 1211 as 



°° 1 f 
S(/3, -04, A) = Si oop = z2j^ 

n=o ' J 



N 



x exp 



JJ dC n z(C n ) 



(2.20) 
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In l|2.20H the notation C = (r, £, a) stands for the loop position r, the loop 
shape £ and the loop species a. When the measure is defined as 

d£=jrfdr[ 2? x , a (0, (2.21) 



A 

simple combinatorics allows one to replace the summation over the N^s 
by a single summation over N = J2 a ■ The loop fugacity depends on 
the distance from the wall as 

z(£) = z a 9(x), (2.22) 

where 6{x) is the unit Heaviside function. The interaction between loops 
arising from the Feynman-Kac formula couples only line elements with the 
same abscissa s 



V{Ci,Cj)= [ dsv(T i + \ ai t i (8)-r j -\ aj t j (8)) 
Jo 



(2.23) 



Thus it is different from the electrostatic potential V e i e ct(£i, £j) between 
uniformly charged wires where any line element of a loop interacts with 
every line element of the other loop, 

V i oct (£i,£j)EE f ds [ da' v ( Ti + \ a M s ) ~ r J - X *M S ')) ■ ( 2 - 24 ) 
Jo Jo 

For a system with quantum statistics, Si oop has the general expression 
written in Ij2.20t where loops L and their fugacities z{L) have more complex 
expressions than in the case of Maxwell-Boltzmann statistics 0]. Quantum 
statistics is taken into account thanks to an extra internal degree of free- 
dom, the number of particles exchanged in the same permutation cycle. 

Equality 12.2011 between the grand partition function of a quantum gas 
of point particles and the grand partition function of a classical system of 
loops with random shapes is the root of an equivalence between both sys- 
tems. As derived in Ref.J], the quantum density p a {x) can be determined 
from the loop density p(£) defined as a grand-canonical average calculated 
with Si 00 p. When exchange effects are neglected, 

p(£) = (^S(r n -r)6($ n -05 an , a \ , (2.25) 

and the relation between particle and loop densities reads 

p a (x) = J V x>a (£) p(C), (2.26) 

where p(£), with C = (r, £, a), does not depend on the projection y of r 
onto the wall plane. 
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2.5 Ideal gas 

In the case of an ideal quantum gas, the grand partition function 12.201 is 
reduced to 

and, by virtue of Q2.25I1 . p ld (C) = z(C) given in 12.221 . The density profiles 
of an ideal gas with fugacities z a 's are given by 112.411 (or equivalently by 
and by using l|2.14j) - l|2.1fijl we retrieve that 



p*(x) = z a 1 - e- 2x + z a O\[^\ . (2.28) 




We stress that, in the vicinity of the wall, the quantum charge fluid 
cannot be handled with as a system made of independent charges, because 
it cannot simultaneously obey the volumic global neutrality H1.13I) and the 
surfacic global neutrality H1.15|) . Indeed, the bulk densities in the ideal gas 
are equal to the z a 's. If the constraint <I2 . 1 If) is arbitrarily enforced upon 
fugacities, though there is no degeneracy with respect to fugacities in the 
purely noninteracting case, the bulk densities satisfy the bulk local neutral- 
ity relation l|1.13fl . However, the global surface charge of the corresponding 
ideal gas does not vanish when species have different masses: according to 
H2.11H and 112. 28j) . it is equal to 

f°° 1 Pir 

a id = dx ^e a p ld {x) = --J-J2e a z a \ a . (2.29) 

•'O ^ V ^ 



3 Diagrammatic representation 

The main lines of the following diagrammatic expansions are analogous to 
the formalisms devised for bulk quantum properties and classical density 
profiles in Refs. A and [5], respectively. 



3.1 Generalized fugacity-expansions 

The equivalence with the classical loop system allows one to use techniques 
originally introduced for classical fluids. For instance, the Mayer diagram- 
matics initially built for point particles can be generalized to objects with 
internal degrees of freedom, such as the species a or the loop shape £. 
A generalized Mayer diagram for the loop density p(C) contains one root 
point C which is not integrated over and N internal points (N = 1, . . . oo) 
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which are integrated over, while each pair of points is linked by at most 
one bond 

f(£ i ,C j )=e- /3e ^ v ^' c ^ -1. (3.1) 
We choose to write the loop-fugacity expansion of the loop density as 



{1 r T N 1 



(3.2) 



The summation is performed over all unlabeled, topologically different, 
connected diagrams G where the root point C is not an articulation point. 
An articulation point is defined by the following property: if it is taken 
out of the diagram, the latter is split into at least two pieces not linked 
together by any bond. (In another diagrammatic representation |4j, which 
is analogous to 13.211 but without the exponential, the root point C may be 
an articulation point.) [J| /] G is the product of the /-bonds in diagram G 
and Sq is the symmetry factor, i.e. the number of permutations of internal 
points C n that do not change this product. 




Figure 1: Mayer diagrammatic representation of the loop density. A white 
point stands for the root loop of the diagram, the coordinates of which are 
not integrated over. Black points are internal loop-points, the coordinates 
of which are integrated over. The /-bond between two loops is represented 
by a line and the weight of a black point associated with a loop Li is the 
loop fugacity z(&). 

At large distances with respect to de Broglie thermal wavelengths, the 
loop potential V(C, £') behaves as the Coulomb potential between the total 
charges of loops, as if they were concentrated at positions r and r'. Because 
of the latter l/|r — r'| interactions, the integrals associated with generic 
diagrams G in H3.2J1 diverge in the thermodynamical limit. 



3.2 Systematic resummations of large-distance 
Coulomb divergences 

The large-distance divergences arising from the long-range of Coulomb po- 
tential can be dealt with by introducing auxiliary bonds and by system- 
atically resumming subclasses of auxiliary diagrams G. The method is a 
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generalization of the procedure introduced by Meeron for bulk quantities 
in a classical Coulomb fluid |ll]. The systematic resummation procedure is 
displayed in AppendixElwhere similarities and differences with the process 
used in Ref.[4| are stressed. 

As in Ref.0], the decomposition into auxiliary bonds relies on the mul- 
tipolar decomposition of the loop interaction. This decomposition allows 
one to exhibit classical screening through the appearance of a screened 
potential <j> arising from the resummation process. It reads 

V (A , Cj ) = V cc (r, - r j ) + V c m (r< , Cj ) + V m c (A , rj ) + V m m (A , C 3 ) . (3 .3) 

where V cc (ri — r^) = u(rj — Tj) is the charge-charge - i.e. monopole- 
monopole - interaction between the total loop charges, and the charge- 
multipole interaction V cm is equal to 



ds 



v (rj 



(3.4) 



with a symmetric definition for V mc . The multipole-multipole interaction 

V mm ig 

V mm (£i, Cj) = V(d,Cj) - v(r t - rj ) - V cm (n, £,) - V mc (A,r,) 

= ds v fa + \i£i(s) - Tj - Xjtjis)) - v (n + \i$i(s) - rj) 
Jo 

-vfa-Tj-XjZfofi+vfa-Tj) . (3.5) 

V c m does coincide with the charge-multipole Coulomb interaction between 
a point charge and a charged wire, whereas y mm is not equal to the 
multipole-multipole Coulomb interaction between two charged wires, be- 
cause it couples only line elements of C and £ with the same abscissa. The 
resummation procedure introduces a screened potential 4> arising from the 
sum of all chains built with the auxiliary bond 
f cc {££) = -f3e a e a ,v{r - r') (see FigEJ. It reads 



0e a e a ,<f)(r,r') = f cc (£,£) 



JV=1 



d£ n z[C n ) 



r%c,£)r c (£,c 2 )...r(c N ,£). (3.6) 



The properties of <j> are studied hereafter. 

When resummed bonds are defined, associated resummed weights and 
"excluded-composition" rules ensure a one-to-one correspondence between 
each class in the partition of auxiliary diagrams G and each resummed 
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0e a e a . <f>{v ) v')=£0 Q £' 




Figure 2: Screened potential </>. 

diagram P*. Contrary to what is done in Ref.[4J, we choose to consider 
nine resummed bonds defined hereafter. The reason is that the choice of 
these nine bonds is associated with renormalized weights which are more 
convenient for dealing with the case of a wall with an electrostatic response 
(see Section 2} than the renormalized weights which appear when only the 
five resummed bonds of Ref.[4j are retained. 
The nine resummed bonds are 



F cc (£,£') =-pe a e a ,<f>(r,r'), 



(3.7a) 



F mc (£,£') = -/3e Q e Q , / ds «£( r + A a £(s),r') - 0(r, r') (3.7b) 
Jo 1 J 

(see FigEl with a symmetric definition for F cm , (1/2) [F cc ] 2 (see FigQJ, 
(1/2) [F cm ] 2 , (1/2) [F mc ] 2 , F cc .F cm , F cc .F mc , and 

F RX (£, £') = | e F"+^'«+F- + F— _ 1 _ pec _ pern _ f mc 

l[F c f- ^ cn f- 2 



__ [_F CC ] 2 — - [i^ cm ] 2 _ - [^ mc ] 2 _ pec pem _ pee pmc I (£^£'y 



As shown in Ref.@], 

F mm (£,£') = -f3e a e al mm (£,£') + W{C, £'), 



(3.7c) 



(3.f 



where </> mm is defined as V mm Q3.5H with 0(r, r') in place of v(r — r'), and 
W is a purely quantum contribution, 



W{C,C) = -[3e a e a . [V(£, £') - V cIcct (£,£')] 



(3.9) 



„1 „1 +oo +OC j ^ 

-(3e a e a > ds ds' [S (s - s') - 1] V" -r — 
x [A a £(s).V r ] 9 [A Q ^'( S ').V r ^' «(r - r') 
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We notice that the argument of the exponential in l|3.7c|) can be written 

F cc + pcm + pmo + pmm = p^ + ^ where 

Select = -j3e a e a , J Q ds f Q ds' <£(r + A a £(s), r' + A a /£'(s')). 



F" 



f mc 



fmc f c 



c O- -O c = c o— o a + c 




Figure 3: Resummed bond F" 



l [jF cc]2 



1 [/CC] 2 



cO=OC = co=OC 



£CC £< 



yc c yc 



/ CC T / c 



*• yc c * 



co — oc 

f cc 



Figure 4: Resummed bond i [F cc ] . (The symmetry factors are not written 
in the figure.) 

Eventually, the Mayer fugacity-expansion i°>.2t of the loop density can 
be rewritten as 



A' 



n ^ 

i=i 



[n 



F 



(3.10) 



The effective screened fugacity arising from the resummation of Coulomb 
ring diagrams (defined in Appendix El and FigEJ is equal to 



where z(L) — 9{x)z a and 



(3.11) 



V c Toud(£) = \J dsj^ ds' [0 - v] (r + A^( S ),r + A Q £(s')) • (3-12) 
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The P* diagrams are defined as the G diagrams in ll.'i.2ll apart from the 
following two differences. First the /-bond is replaced by the nine .F-bonds. 
Second, P* diagrams obey an "excluded-composition" rule associated with 
the fact that all points have not the same weight w(C) (in order to avoid 
double-counting) , 

!z sc (£) — z(C) if C is involved only in a 
product F ac (£i,£)F cb (£,£j) (3.13) 
z sc (C) otherwise. 

In H3.13H superscripts a and b stand either for c or m, and the points Li 
and Cj may coincide. 




Figure 5: The sum of ring diagrams attached to a white point, I T = 
— /3e Q Vpf oud . 

The vicinity of the wall replaces the bulk exponential screening by an 
integrable algebraic screening in the directions parallel to the wall (see next 
section), and all diagrams with resummed bonds are finite in the thermo- 
dynamical limit, when integrations are performed first over loop shapes 
and then over loop positions. The reasons are the following ones. First, 
in F cc l|3.7afl the resummed charge-charge interaction between the total 
loop charges is proportional to the screened potential <j>(r, r') which obeys 
an inhomogeneous Debye equation where the screening length depends on 
the distance from the wall, as discussed in next section. The translational 
invariance along the wall ensures that 4>(r,r') = 4>(x,x',y) where y is the 
projection of r — r' onto the wall. As shown in an analogous classical situ- 
ation [12j . where the inhomogeneity in the zero-coupling limit arises from 
the electrostatic response of the wall, when |y| goes to infinity while the 
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values of x and x' are kept fixed </>(r,r') decays as 



0(r,r') r 1^1. (3.14) 

|y|->+oo |y|° 

Far away from the wall the x-dependent screening length tends to a non- 
zero value, so that f(x, x') decays exponentially fast to zero at large x or x' . 
Therefore the leading tail of F cc has the structure 
g cc (x, x')/\y\ 3 and F cc is integrable. In F nlc l|3.7bj) the resummed 
multipole-charge interaction is also integrable. Indeed, since the Brown- 
ian path measure T> x ,a(£) ensures that all moments of £ are finite, and the 
leading tail in F mc , which is not canceled by the integration over the loop 
shape £, behaves as 



.1 +oo 



, S10 d q g cc (x,x') 1 



d q g cc (x, x')/dx q is an exponentially vanishing function of x and x' far away 
from the wall, while J Q ds T> Xja (£) [t; x (s)] q is a gaussianly vanishing function 
of x over the scale A Q . As a consequence, the bond F cm is also integrable 
in x and x' over the scales A a and X a > , respectively. By virtue of Il3.7c|l the 
tail of F RT (£,£') at large distances is given by the tail of the resummed 
multipole-multipole interaction in F m m 113.811 and can be decomposed into 
a "diffraction" term arising from </> mm and a purely quantum contribution 
described by W defined in 13. 9J1 . The corresponding leading tails are 

r l ,.\ +oo +oo 



.! 1 t ,„m<j t> ,jsW d q+q 'g cc (x,x') ^ 1 

q=lq' = l ' 



r r J—Z _ _ i i 



(3.16) 



and 



-><..<„< / -A- / ' ' ' ' (3 ' 17) 



o Jo 



The tail H3.16I) is integrable for the same reason as the tail H3.15I) of F m c . 
When |y| goes to infinity while the values of x and x' are kept fixed, the 
function l|3.17JI decays as |y| times a function which, after integration 
over the Brownian measures T> x . a (£) and I? x / )7 (£ ), depends on x and x' 
and converges gaussianly fast to zero over the scale A Q and A Q ' respectively. 
Eventually, thanks to resummations, every diagram in H3.10I) is well defined 
in the thermodynamical limit, and we will consider this limit from now on. 
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4 Screened potential 

4.1 Debye equation for an inhomogeneous fluid 

The screened potential 4> which is defined as the sum of chains \?>M is the 
solution of the integral equation 

<j>(r,r')=v(r ) r')-±- j dv"K\x")v{v,v")^{v",v'). (4.1) 

In 114,111 the positive function ~k 2 (x) reads 

k\x) = 4^/3^4 f V x , a {£)z{C) 

a •* 

= 6{x) 4vr/3^ e 2 a z a [l - e" 2a;2 / A =] , (4.2) 

a 

where the second equality arises from H2.14Jl - H2.16ll . Far away from the 
wall k 2 (x) tends towards its bulk value 

47r/3^ e 2z Q . (4.3) 



k 2 



By using Poisson equation satisfied by Coulomb potential v, 

A r v(r, r') = -47r5(r - r'), (4.4) 

the integral equation i|4.1H is shown to be equivalent to a set of partial 
derivative equations, the explicit expressions of which depend on the signs 
of x and x' . For x' > 0, 0(r,r') is a solution of 

A r 0(r, r') - k 2 (x) <p(r, r') = -47n5(r - r') for x > (4.5a) 

and 

A r 0(r,r')=O for x < 0. (4.5b) 

Moreover, the diagrammatic definition i|3.6[l of 4> implies that 4> obeys the 
same boundary conditions as the electrostatic potential v(r,r'), 

c^(r,r') 



(r,r') and 



Ox 



are continuous at x = (4.6a) 



lim <t>(r,r') = 0. (4.6b) 

r| — >oo 

The invariance of the system in directions parallel to the interface im- 
plies that the Fourier transform of 4> along these directions obeys a one- 
dimensional differential equation with respect to x. We introduce the di- 
mensionless coordinates x = kx, y = «y and r = nr and the dimensionless 
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screened potential </> defined by 

$(x,x',y) = ~c/>(x,x',y). (4.7) 

K 

The dimensionless Fourier transform of <j> along directions parallel to the 
wall reads 

^(x,x',q) = J dye-^{x,x',y). (4.8) 
For x! > it is a solution of the differential equations 



dx 2 



and 



(l + q 2 ) -U{x)\^(x,x',q) = -Att5{x - x') for x > (4.9a) 
-q 2 |0(x,J',q) = for£<0, (4.9b) 



with the potential function 



U(x) = - E " e ^ e " 2 . (4.10) 

4.2 Series representation of the exact screened poten- 
tial 

It is well known (see e.g. Ref. ^21) that the solution of the one-dimensional 
differential equation i!4,9a|) can be formally written in terms of the solutions 
h of the associated "homogeneous" equation (where the Dirac distribution 
on the r.h.s. is replaced by zero and) which is valid for — oo < x < +oo, 

{^-( 1 + q 2 )-^)}M^;q 2 ) = 0. (4.11) 

More precisely, a solution (j> of l|4.9ajl can be decomposed into the following 
sum: a linear combination of two independent solutions h + and h~ plus a 
particular solution </> S ing which is singular when x — x' and is calculated in 
terms of h + and h~ by the so-called Wronskian method. In the following, 
h + (h~) is chosen to be a solution which vanishes (diverges) when x tends 
to +oo. Then s i ng reads 

sing (£,£',q 2 ) =- ^— /t~(inf(x,.T , );q 2 )^ + (sup(x,x , );q 2 ), (4.12) 
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where the Wronskian W^q 2 ) = h (dh + /dx) — h + (dh /dx) is independent 
of x. Since the differential operator in (I4,9aj) is self-adjoint, the real poten- 
tial 4* has the symmetry (j)(x, x',q) = (f)(x',x, q) when x > and x' > 0. 



Then the boundary condition i|4,6b(l obeyed by <ft enforces that 

4(x, x',d) = s ing(i, x, q 2 ) + Z(\q\) h+(x; q 2 ) q 2 ), (4.13) 

where Z(|q|) is determined by the continuity relations Ij4.fi ajl . 

The particular solution h + (h~) of the "homogeneous" equation (14 . X If) 
can be written as the solution exp[— x\J\ + q 2 ] (exp[x^l + q 2 ]) of (14. 1 If) 
when U = 0, times 1 + H + (1 + H~) where the function H + (H~), which 
describes boundary effects, is chosen to vanish at x = 0, 



h ± {x;q 2 ) = eT*V I +^[l + H ± (x,q 2 )}. (4.14) 

Then 

6 ■ (x x' a 2 ) - i?!_ e -|5-5'lv / I+? 

x [1 + H- (mi(x,x');q 2 )][l + H+(sup{x,x');q 2 )}. (4.15) 
The Wronskian W{q 2 ) proves to take the simple form 

(4.16) 



w(q 2 ) = -2VYT^ + d -^§^l 



.7—0 



Indeed, H + (x;q 2 ) and H~(x;q 2 ) vanish at x = 0, and so does dH~/dx, 
as can be checked on the formal solution given in next paragraph. 

As shown in Ref.|2J, H + can be written as the formal alternative series 

H+(x;q 2 ) = -T+[l}(x;q 2 )+T + [T+[l]] (x;q 2 ) (4.17) 
where the operator T + acting on a function / reads 

pX Z' + OO 

T+[f](x;q 2 ) = dve^y 1 ^ dte~ 2t ^+^ U{t)f{t). (4.18) 

JO Jv 

The solution H~ can be similarly written in terms of another formal series 
ff-(5?;q 2 )=T-[iP;q 2 )+T- [T"[l]] (x;q 2 ) + -.- (4.19) 



with 



^[.f]{x:q 2 )= I dve~ 2v ^+^ I dt e 2t ^+^ U (t) f (t) . (4.20) 
Jo Jo 
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By combining I|4.13|l - lj4.2()|l . we get a representation of the screened poten- 
tial 4> in terms of formal series. 

These series are bounded by geometric series of kX (with A = sup a {A a }), 
because the Gaussian function U is integrable for all x's. More precisely, a 
straightforward calculation shows that 



T+[1P) 



'\ a y/T+^) 2 /2 



; 2 2 yiw Erfc f vjx + Wi + q : 



A a 



V2 



-Erfc 



V2 



(4.21) 



where A a = n\ a , and this expression obeys the inequality 

|T + [1](2)| < kX. 



(4.22) 



4.3 KA-expansion of the screened potential 

In Sectional we will restrict our explicit calculations to a low-degeneracy 
and weak-coupling regime (see ljl.41) and l|1.7ll V In this regime ft A is also 
negligible with respect to 1 : the de Broglie thermal wavelengths A Q 's are 
small compared with the typical screening length k . 

Let 4>^ be the solution of equation (l4.Qa.fl when the profile k 2 (x) is 
replaced by its bulk limit k 2 . The^aboye formal series for H + and H~ 
provide a systematic expansion of <ft — <ft^ in terms of the ratios of the 
length scales A Q 's, over which k 2 (x) varies, and the length scale K , which 
is the limit of k~ 1 (x) when x goes to infinity. The expansions of H + and 
H~ are absolutely and uniformly convergent with respect to the variable x 
for values of kX smaller than some finite value, because they are bounded 
by geometric series. As a result, 

x', y) = 0<°> (x, x',y)+ 0(kX), (4.23) 

where the leading-order term <//°) is equal to (r, r') (see H4.7J1 1. By 

definition, </>(°)(r, r') is the solution of 

A r ^°> (r, r') - 9{x) n 2 0<°> (r, r') = -47r<5(r - r'), (4.24) 

with the same boundary conditions as </>(r,r'). According to 114. 1311 . we get 
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Vl + q^ { y/1 + q 2 + q J 

(4.25) 

Moreover, by definition of the screened potential, the correction O(kX) 
in <I4. 23|) is independent of the root-point species a and a' . As shown in 
Section 15.31 the latter property implies that the contribution to density 
profiles from this 0(KA)-correction in </> is canceled at first order in kX by 
the neutrality constraint 12. lit on fugacities. It is the reason why we do 
not write the explicit expression of the term O(kX) in H4.23JI . 



5 Weak-coupling expansions 

5.1 Subregime for parameters e and (A/a) 3 

In the weak-coupling and low-degeneracy regime ljl.71) and H1.4J1 . only a 
finite number of resummed Mayer diagrams in the representation H3.10I) of 
the loop-density p{£) contribute at lowest orders in the classical coupling 
parameter e = (l/2)n(3e 2 . Moreover diagram contributions may be also 
expanded in powers of kX. 

Indeed, kX is a function of the two independent parameters (A/a) 3 and 
e: by virtue of 1)1. fijl . 

kAoc^£ 1/3 , (5.1) 

and kX <C 1 in the considered regime lll,7|) and 11.411 . Then the screened 
potential (f>, which is involved both in screened fugacities and in resummed 
bonds, can be expanded in powers of kA, as well as the functions resulting 
from integrations over Brownian paths. As shown by the scaling analysis 
performed in next Sections 15.21 and l5~3l the diagrammatic representation 
of the particle density derived from 12.2(it and ffi.lOjl provides a systematic 
expansion of p a {x) in powers of parameters e and kX, where exchange 
effects are neglected. 

More precisely, we will show that, because of quantum dynamics, the 
first coupling corrections to the ideal-gas particle density is a sum of two 
terms of order e and kX respectively, and the next contributions are of 
order 

e\ e 2 |ln( K A)|, e.nX, (kX) 2 , f^j Q w (-^) , (5-2) 



with 



A ^oc^! and ^oc^ (5.3) 
e X kX 
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Q w {t) is expected to vanish at t — by analogy with the bulk function Q 
defined in Ref.|14|, On the other hand, exchange corrections, which have 
been neglected from the start, are of order 



A 



where E w (t) is expected to vanish as t goes to zero, for the same reason as 
Qw(t). 

As a consequence, we shall be allowed to retain only the corrections 
linear in e and kX in p a {x), if these terms are larger than the exchange 
corrections of order ll-').4ll and the coupling corrections of order i!5 .3fl . We 
consider the subregime where 

A x 3 



£ < y-j « £, (5.5) 

In 15. ot e 2 < (A/a) 3 means that the ratio of e 2 and (A/a) 3 is either kept 
fixed or tends to zero when both e and A/a vanish. In other words, /3e 2 /X 
is kept fixed - i.e. the temperature remains fixed - or vanishes - i.e. the 
temperature goes to infinity. In this subregime the condition 

A\ 3 „ / (3e 2 



F ( -'— I « e, (5.6) 
with F = Q w or E w , is met, as well as the condition 

F (-(3e 2 /X) < kX. (5.7) 



A^ 3 



Similarly, inequalities e 2 <§; nX, e 2 |ln(KA)| e, e 2 |ln(KA)| <C kX, and 
(/tA) 2 <C £, are also satisfied. Eventually, in the subregime 115. 5J1 . which 
can be rewritten as e 3 < (kA) 3 -C £ 2 by virtue of llo.lll . we may retain 
only contributions of order e and kX, and the neglected terms are of order 
O (?7 2 ), where rj 2 is a generic notation for the terms in H5.2J1 and 115.411 . 

<D(rj 2 ) = C((|Ol, (El)- (5.8) 



5.2 Screened loop fugacity 

In this section we calculate the «;A-expansion of the screened fugacities 
ffi.lli . As shown in Ref. [JJ, the free energy e 2 Vpf oud is associated with the 
"geometric" repulsion from the wall due to the deformation of the screening 
cloud surrounding every charge near a boundary. According to l|3.12fl and 
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the fact that (</> — v) and its derivative are continuous (because 4> an d v 
have the same singularity when r = r'), 

V c Toud(^) = Kioud(r) + -O(kX). (5.9) 

K 

H5.9I) is the Taylor expansion of V^ c oud around its classical value 

VZ ud (r) = ^[cf>-v](r,r). (5.10) 

(We stress that V r c s 1 Q Ud (r) has no singularity in the range < x so that no 
classical spurious singularity is introduced by the expansion Il5.9|) .l The 
KA-expansion (|4.23|) of the screened potential (j> leads to 



Kfoud(r)=^(r) + ^(«A), 



where at first order 

„2 T/ sc(0) 



-[3e a V c 



2 r 



loudV r ) — 2 



In H5.12J1 we have used the definition 



(5.11) 



(r,r)=e a [l-£(/«c)]. (5.12) 



(5.13) 



By virtue of Il4.25jl . 

L(u) 



oo 

1 

-2u 



dt 



ill 

2u u 2 2u 3 



-K 2 (2u), 



(5.14) 



where K2{2u) is a Bessel function, which decays proportionally to 
exp(— 2u)/^/u at large u. L(u) is a continuous positive decreasing func- 
tion for u > 0. Therefore L(u) is bounded by L(0) = 1/3 and, since e a is 
small, we can expand the exponential function in the definition H3.11|) of 
screened loop fugacities. We get 



z sc (C) = 9{x) z a \ 1 + e a [1 - L(kx)] + Q(s . k\) 



(5.15) 



Since L(u) vanishes exponentially fast when u goes to infinity, z sc {C) tends 
to z a {1 + £ a } far away from the wall. 
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5.3 Diagram contributions at leading orders 
5.3.1 Bond F cc 

The integral associated with the diagram with one F c c -bond reads 

J dC'z sc (C')F cc (C, £') = -/3e a ^ e 7 J dr'<f>(r, r') 3^V), (5.16) 



where the screened fugacity 2^ c (x) is defined as 

z%(x) = J V x . a {i)z sc {C). 



(5.17) 



By using the value II5.15|I of the screened loop fugacity and the integral 
H2.16fl of the path measure, we get 



l-e 



1 + e a [1 - L(kx)] \ + z a O{e.n\). (5.18) 



By virtue of i!4,23JI . the leading-order terms in i|5,16|) arise from the struc- 
ture 



-fie a ^2 ' 



dr' 



x'>0 



^°)(r,r') + 0(^)1 \l-e- 2x ' 2 / x "] [l + 7 (e) 



(5.19) 

where 7 (e) is a term of order e which depends on the species 7. (0 7 (e) 
is the screened self-energy term in ~z^(x).) 

The a priori leading-order term in Q5.19H is obtained by retaining only 
the two constants 1 in brackets and <^°) (r, r') , This term is independent 
of the de Broglie wavelengths, because (/>(°)(r, r') is purely classical and 
involves only the length scale 1/k d . By virtue of H4.25J1 . 

/ dr' (r, r') = — / dx' { e — 12— 2'| + e -(S+2')l ( 52Q ) 

Jx'>0 K JO ' 

As a consequence, the leading term in H5.19J1 is both O(e ) and O ((kA)°), 
namely, before summation over 7, 

J dr' J 2V, 7 (0 z sc (£')F cc (£, C) = Oil). (5.21) 

However, after summation over species, the 0(1) contribution in H5,lfi|l . 
where the quantity J2 7 e 7 z 7 is factorized, is exactly equal to zero because 
of the neutrality condition H2.11I) imposed on fugacities. Moreover, with- 
out condition H2.11I) . we would have to consider an infinite number of 
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diagrams at leading order, because the addition of a "star" subdiagram 
nr=2 [/ dCiZ sc (Ci)F cc (C , , with an arbitrary number n, to the di- 
agram J d£'z sc (£')F cc (£,£') would also yield a contribution of leading 
order 0(1). Nevertheless, we notice that, after summation over all dia- 
grams, the final expansion of the density must be independent of whether 
the condition on fugacities is fulfilled or not, because of the degeneracy 
among fugacities discussed in Section l2~2l 

In fact, the diagram with one bond F cc contributes at orders e and kX, 



The contribution of order 0(e) comes from </> (0) (r,r') X 7 (e) in {5T9t . 
and it has been calculated in Ref. jH] ■ The contribution of order kX arising 
from the product of constants 1 in ljfi.1911 times the term of order kX in the 
expansion of 4>(r,r') is canceled by the neutrality condition H2.11JI , (This 
is also true for the whole f«A-expansion of 4>{r, r'), because all terms which 
depend on a species 7 only through the product e 1 z 1 are canceled when the 
summation over 7 is performed.) Therefore, the contribution of order kX 
arises only from <jr- ' (r, r') exp[— 2cc' /A?,] in Il5.19|) , The factor kX is yielded 
by the ^'-integration of the Gaussian term which arises from the integrated 
quantum measure l|2.16JI . Indeed, for any bounded and integrable function 
/ which decays less quickly than exp[—2x 2 /X 2 t ] at large x 



ndx e- 2x ' 2/x «f(K,x) ~ K X a x /(0) / dte~ 2t2 (5.23) 



(The Gaussian factor in the integral on the l.h.s. of l|5.23|) makes this 
integral convergent over the scale A 7 , and not k _1 as it would be the case if 
this factor were not here.) As shown in Appendix El the same mechanism 
also operates for the ^'-integration of odd moments of a Brownian path 
£' (involved in diagrams with F cm -bonds), because these moments decay 
gaussianly fast to zero at large x' over the scales A Q 's (see for instance 
l|2.18|l ). contrary to even moments which tend to their non-zero bulk values 
away from the wall. 

5.3.2 Other resummed bonds 

As shown in Appendix El by combining previous arguments, Taylor ex- 
pansions around classical expressions and construction rule Q3.13J1 about 
fugacities, where z sc (£') — z{£') — z Q 0(e), a simple scaling analysis shows 



that diagrams with one bond, ,F cm , F mc , [F cc ] 2 /2, F cc F cm , [F cm ] 2 /2, 




(5.22) 
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F cc F mc , or [F mc ] 2 /2 are of orders equal to either {kX) 2 , e ■ kX, s 2 , s{kX) 2 , 
e 2 ■ k\, or e 2 ■ (kX) 2 . The order of the contribution from the bond F RT is 
determined by analogy with its order in the bulk case. 

5.3.3 Global results 

As a result, the complete scaling analysis shows that, at first order in 
e and kX, the loop density comes from l|3.1f)|) where only one diagram, 
namely the diagram with a single F cc -bond, is retained in the argument of 
the exponential. Moreover, the expression of this first-order contribution 
involves only the lowest-order term <fr(°> (r, r') in the KA-expansion H4.23JI of 
the screened potential. The loop density reads 

p(C)=z sc (C) 

xexp j-/?e Q ][> 7 y drV>(°>(r,r') (^J © x /, 7 (0 2 SC (£')) 

+ (5.24) 



F cc 

p(C) = z sc (£) exp { CO • z sc (C) + 0( V 2 ) 



Figure 6: Diagrammatic expression of the loop density at first order in e 
and kX. 

Since the argument of the exponential in ljo.2411 proves to be a bounded 
function of order e and kA, the loop density is given at first order in e and kA 
by linearizing the exponential in H5.24J1 . According to H2.26J1 . the particle 
density profile is obtained by performing the path integration J T> x ^(£) 
with the result 

Pa (x) = z s a c (x)^l-/3e a ^2e 7 J dv' z^(x') 4> (0) (r, r') + 0{if) | , (5.25) 

where z^ix) ls gi ven m l|5.18|l and rj 2 is defined in Ij5.8ll . 
5.4 Electrostatic potential 

In this section we show how the integral involved in the expression H5.25J1 is 
related to the potential $(x), which is defined as the difference between the 
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electrostatic potential created by the fluid and the electrostatic potential 
V R in the particle reservoir (located in the bulk). This potential obeys 
Poisson equation 

— (a:) = -4^}Je a( o a (x) (5.26) 

a 

with the boundary conditions : $ and d^/dx tend to when x goes to +00. 
The condition about the derivative of <I> arises from the absence of any 
net electrostatic field in the bulk for a Coulomb fluid at equilibrium. The 
condition about $ fixes the electrostatic potential reference. The definition 
of leads to the integral representation 

/+00 
dx'(x' - x)^e a p a {x'). (5.27) 
a 

As in the classical case 0, the structure H5.25JI of density profiles can 
be rewritten in the form 



p a [x) = z%(x) l-Pe a G[x) 

with 



(5.28) 



G(x) = J2 e ^ I dr'^ c (x')^°Hr,r l ). (5.29) 
7 J 

As checked in Section l5~5l at leading order G{x) is a function of kx. There- 
fore, according to II5.18|I and (I5.23|) . the contribution from 
— P e a \^ai x ) ~ z a\ G{x) to the integral in H5.27J1 is a correction of rel- 
ative orders O(kX) and 0{e) with respect to the leading contribution from 



a \ a / 



(5.30) 



Therefore, the argument of Section 5.4 in Ref. holds. It reads as fol- 
lows. According to its definition and the partial derivative equation 114.2411 
satisfied by tfi^ , G[x) obeys the differential equation 



d 2 G(x) 
dx 2 



-Ait^e a !^{x)+K 2 G{x). (5.31) 



Since dG(x)/dx is finite and decays to zero when x goes to infinity, combi- 
nation of i|5.27|l . (I5.30|) . and H5.31J1 implies that the electrostatic potential 
is merely 

G(x)- lim G(x) ' 



x— »+oo 



+ — 0(e, K \), (5.32) 
pe 
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where 0(e, kX) denotes a sum of terms of order e and kX respectively. 

Eventually, by virtue of l|5.18|l and Q5.32I1 , the density profile 115, 28|) can 
be rewritten in terms of the bulk density and the electrostatic potential 
drop created by the fluid as 



p a {x) 



1 



{l - e a L( K x) - 0e a $(x)} + p B a O{rf), (5.33) 



where the bulk density is given by 

Pa 



lim p a {x) = z a < 1 + e a — (3e a lim G(x) + 0(r\ ) 

x — >+oo 



(5.34) 



Since the bulk density p B coincides with z a at leading order, k d defined in 
il .ol is also equal to k at leading order 



1 + O(e) 



and e n = e 



1 + O(e) 



(5.35) 



This will enable us to consider the Debye screening length as the refer- 
ence length scale for classical effects when writing the final results in next 
section. 



5.5 Decomposition into classical and quantum contri- 
butions 

In the density profile H5.33J1 . the term involving L is purely classical, whereas 
the electrostatic potential can be split into a classical contribution $ cl ( e ) of 
order O (e/(/?e)) and a quantum contribution $i u ( kA ) of order O («A/(/3e)), 

= $q u (« A ) ( x ) + $ cl ( £ ) (x) + 4-0 (rf) . (5.36) 

pe 

Accordin g to l5~53f, $ cl ^(x) = G cl( - E \x) - lim w ^ +00 G cl ^ (x) with, by 
virtue of II5.29J1 , 115,1 8J1 , and of the neutrality constraint (12 . 1 lfl upon fugac- 
ities, 

G cHe \x) = ^e 7 f dr'9(x')z^ [l -L(kx')] </> (0) (r,r'), (5.37) 

7 J 

whereas $i u ( kA ) is the leading term in the expansion of 

- e 7 / dr' 6(x') z 7 e- 2x ' 2/x2 -< (o) (r, r'). (5.38) 

7 ■* 
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Since J dytj)^ (x, x',y) is a function of kx' which decays exponentially fast 
when kx' goes to infinity, according to l|5.23ll , 



j^dx' (^-Y^e^e-^' 2 /^ Jdy^(x,x' = 0,y). 



Jo V '7 

(5.39) 

The classical part $ cl ( £ ) in the electrostatic potential has already been 
calculated in Ref.Q] with the result 

$ cl(e) (x) = -AM(k d x), (5.40) 

where the function M is 

r°° e ~ 2tu - 2te~ u 1 

"<•>-/ dt i-w ( » + v ^i)- (5 - 41) 



and the constant A reads 



T e 3 n B 

A=^4^==- (5.42) 



In H5.4()|l . the argument of M has been written k d x in place of kx, by virtue 
of (fo~35| . 

The purely quantum part of the electrostatic potential drop is derived 
from II5.39|I . According to the expression of / dy(f>^{x,x' , y) already used 
in 

$ qu( re A)( x ) = _fr,g e -^ ; (5.43) 
where the constant B depends only on the fluid composition 

„ (5 . 44) 

We notice that the total electrostatic potential drop between the wall and 
the fluid bulk (set as the reference of electrostatic potentials) is equal at 
leading order to 

" - 1 \ - 



$(0) = --(ki3-l---^J-ftB. (5.45) 

$(0) includes both classical and quantum corrections. 

The structure of quantum particle densities derived from H5.33JI , H5.40JI 
and II5.43J1 is summarized in Ijl.lOJI . In the fugacity expansion of the bulk 
density the explicit first-order correction to the value z a in an ideal 
gas is derived from H5.34JI and the explicit value of G(x). Since 
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tends towards zero when x goes to infinity, the relation between p^ and z a 
does not include quantum contributions proportional to ti and reads 



e 3 « B 

7 7^7 



e aPa 



(5.46) 

Indeed, in the bulk, the spherical symmetry enforces quantum-dynamical 
coupling effects to be at least of order (kA) 2 , proportional to ti 2 , whereas 
exchange effects are at least of order (A/a) 3 , proportional to ti 3 . The ex- 
pression IJ5.46I) is in agreement with the result (5.28) of Ref.pj] calculated 
directly in the bulk. It does satisfy the electroneutrality condition ll.l^jl . 

The expression ll.lOjl of the quantum density profile at first order in e 
and k\ can be rewritten in terms of the densities p^'s at first order in e in 
the corresponding classical system, 



1 - e a L(K D x) - /3e Q $ cl(£) (x) 



(5.47) 



where L is defined in H5.14J1 . At this leading order the classical density 
profile does not involve the short-range repulsion that must be introduced 
in order to prevent the collapse of the system in the limit where h tends to 
zero (see e.g. Ref.pQ). Indeed, Il5.47|) is obtained in a subregime where the 
range a of the short-distance repulsion is such that e 2 < (c/a) 3 <§; e. We 
get 

p a (x) = \l - e-^A-l L^ e \x)-p B a (3e a ^ KX Hx)\+p B a O(r 1 2 ). (5.48) 



The latter expression displays two quantum effects. We stress that the 
effect linked to the vanishing of wave-functions has an essential singularity 
in h. The quantum contribution linear in H in the electrostatic potential is 
allowed by the breakdown of spherical symmetry, whereas bulk quantum 
effects in the physical regime of interest appear only at order ti 2 jjj], as 
already mentioned. 



6 Generic properties 
6.1 Density profiles 

The structure ll.lfljl of density profiles is ruled by the competition between 
three effects. The purely quantum contribution 1 — e~ 2x / A ° arises from 
the vanishing of wave-functions inside the wall. In the physical regime of 
interest this effect is the same one as in an ideal gas (see l|1.2fll. up to 
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amplitude corrections arising from Coulomb coupling. The second term, 
involving the function L, describes the geometric repulsion due to the de- 
formation of screening clouds near a wall pp. Indeed, a charge and its 
surrounding screening cloud are more stable in a spherical geometry than 
in the dissymmetric configurations enforced by the presence of a wall. This 
effect is purely classical at the order of the present calculation. The term 
e a $(x) describes the interaction between a particle with a charge e a and 
the electrostatic potential drop, created by the fluid itself, with respect to 
the bulk (set as the reference of electrostatic potentials) . This contribution 
contains both quantum and classical effects. 

In the very vicinity of the wall, for x < A <C £, D (with A = sup a {A Q }), 
densities are mostly ruled by the quantum effect of the cancellation of 
wave-functions inside the wall, 



Pa(x) ~ Pi l-e- 2 * 2 / A » \ 1 - ^K D Pe 2 a - /3e Q $(0) }, (6.1) 



1 



G 



where $(0) is given in H5.45I) . The heavier a particle species is, the steeper 
the vanishing of its density occurs, since dynamical quantum effects are 
less important for heavy particles. Densities and their first derivatives 
are continuous on the wall, as expected, since densities involve the squared 
moduli of wavefunctions and the latter ones are continuous at the boundary 
of a wall with a possible step variation in their first derivatives. 

At distances from the wall large compared with the quantum de Broglie 
wavelengths, x 3> A, densities vary over the classical Debye screening length 



-n D f3e 2 a L(n D x) + (3e a 



AM(n D x) + hBe~ 



(6.2) 

where A and B are given in 115. 42|) and H5.44H respectively. In this region, 
density profiles are determined by the interplay between the effect of the 
classical geometric repulsion described by L(k d x) and the effect of the elec- 
trostatic potential, with both classical and quantum origins. L{k d x) and 
M(k d x) vanish exponentially fast over the scales £ D /2 and £ D respectively. 

As a consequence, at distances from the wall large compared with the 
classical Debye screening length £ D , the contribution from the electrostatic 
potential dominates in the density profiles 



Pcc{x) - Pa 



where 



Pa *»« 

A 



-/3e Q $ as e" 



ln3+ —= - 2 
V3 



KB. 



(6.3) 



(6.4) 
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Equation 116.311 is valid in the generic case where <& as ^ 0. According to 
<5.42t and Q5.44I1 . <£> as is likely to vanish only in a two-component plasma 
where charges are opposite, e_ = — e + (A = 0), and where both species 
have the same mass (B = 0). When <£> as ^ 0, if the charge e a has a sign 
opposite to that of <& as , Pa{x) > p% at sufficiently large distances x, as 
shown by 116.311 , 

6.2 Profile of the total particle density 

At the order of calculations, the wall is repulsive for the global particle 
density everywhere in the Coulomb fluid, as in the classical case pQ, 

X>aC0<XX (6-5) 

a a 

Indeed, according to the structure of densities, since p a {x) and 

1 — exp [— 2x 2 /A^] are positive, the second factor on the r.h.s. of is 
also positive in the considered regime of small parameters. Therefore the 
effect of the vanishing of wavefunctions near the wall is to lower the density 
p a {x): H1.10|l implies that at any distance x from the wall 

a a 

The contribution from the electrostatic potential drop to the bound 
in ijfi.fiH vanishes because of the bulk electroneutrality l|1.13|l , as in the case 
of Pa^ (x) . Thus the bound involves only the sum of the contributions 
from classical screened self-energies. The corresponding geometric repul- 
sion from the wall tends to reduce the density of each species with respect 
to its bulk value, and we get ljfi.511 . 

6.3 Charge density profile 

Even if the Coulomb fluid remains globally neutral, when species have 
different masses, the local charge density ^ Q e a p a (x) is non zero, as well as 
the associated electrostatic potential drop $(x). The property holds even in 
the case of a charge-symmetric two-component plasma where the classical 
charge density J2 a e aPa( x ) vanishes for symmetry reasons (because the two 
species have opposite charges). (The latter cancellation can be checked at 
first order in e where, by virtue of 115. 4711 . ^ Q e a pa^ is proportional to 
S 7 e 7/°7-) The charge density profile is organized in various layers with 
opposite signs which depend on the composition of the fluid. 

In the very vicinity of the wall (x < A), when species have different 
masses, the charge-density profile exhibits a zeroth-order effect arising from 



(3e 2 a L(k d x) - /3e a $(x) 



(6.6) 
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the cancellations of the various wave-functions over different scales. Indeed, 
according to llfi.lll and l|1.13ll , 

5> Q p Q (z) b< ~ -^Te^e" 2 ^ +epO(e,n\), (6.7) 

where p is the typical particle density. When all masses are equal, 
J2 a e aPa(x) is only of order epO (e, kX) in this region, by virtue of the 
bulk local charge neutrality. 

At distances from the wall large with respect to the quantum lengths, 
x ^> A, the charge density is an effect of order epO (e, kX) which is ruled 
by the competition between the geometric repulsion from the wall and the 
electrostatic potential drop (see H6.2J1 1, According to H5.48J1 and Q5.43I1 . 

J2^ Pa (x) ~ J2e a p c ^ ) (x) + hB^e-^, 

t— 1 A<z<f D 47T 

a a 

where, by virtue of H5.47jl and H5.40I) . 

^2e a p& e \x)=-±K D (j2 e iP°) [L{k d x) -M{ Kd x)] . (6.8) 

a \ 7 / 

In the case of a charge-symmetric two-component plasma ^ a e- a p a {x) is 
purely quantum at distances x X. 

Eventually, the large-distance behavior of ^ a e. a p a (x) is merely ruled 
by the electrostatic potential drop 3>(x), which includes both quantum and 
classical effects, 



y^^gPajx) 



z»Cd 47T 



<$> c 5 e) - m 



(6.9) 



where is given in H6.4J1 ■ As it is the case for the electrostatic potential 

drop and for particle densities, the charge density includes quantum effects, 
linear in h, which exist far away from the wall over a few Debye screening 
lengths. 

6.4 Global charge 

The wall that we consider does not carry any external surface charge and 
the global surface electroneutrality ljl.1511 of the Coulomb fluid is fulfilled, 
as checked in the present section. At leading order O (e, kX) the global 
surface charge a can be decomposed into leading classical and quantum 
contributions: a — a cX ^ + a qu ^ x \ where 

a cl ^ = / dxJ2e a p c M(x), (6.10) 
Jo 
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and cr qu ( KA ) is the term of order kA in 

r + 00 



dx e n 



a [ Pa ( X ) - pf S Hx)] . (6.11) 



We recall that, as mentioned after Q5.47I1 . the hard-core that must be intro- 
duced in order to prevent the classical collapse between opposite charges 
does not appear at order e in p c a (x). 

Since the classical densities already obey H1.15J1 . which is enforced by 
macroscopic electrostatics, a cX ^ and cr qu ^ A ^ must vanish separately. As 
checked in Ref.JT], the classical contribution a cl ^ of order e arising from 

E e aPa ( ^ (x) does vanish. 

The quantum term cr qu ( KA ) of order kA arises only from the two quan- 
tum terms -p B a exp [ -2x 2 / A^] and -pi (3e a <S> qu ( KX \x) in p a {x) - pt [E \x). 
Indeed, according to Q5.48I1 . 

P a (x) - p$ s \x) = -p B a e- 2 * 2 / x « - p B a 0e a ^ KX \x) (6.12) 
- { [pf e) (x) P B a ] p B a Pe a <S>^(x)} e" 2 ^. 

The term in curly brackets is a function of order p„C(e, kX) which decays 
exponentially fast over the scale K ; by virtue of Q5.23I1 . after multipli- 
cation by exp [— 2x 2 /A 2 ] this term gives a contribution of relative order 
O (e . kX, (kA) 2 ) to cr qu ( KA ' defined in ilfTTT)) . Therefore 

cr qu(«A) = CT qu(«A) + ^u( K X) ^ 



and 



qu(ftA) _ 



/•+00 

/ ^^e 2 p-/3$ qu ( KA )(x), (6.14) 
Jo 



where $q u ( KA )(iz;) is given in H5.43I) . 

The ^.-contribution cr^ u ^ A ' 1 to the global charge a is canceled by the 
contribution <j^ kX \ Therefore, the quantum term /3e ct $ qu ^ A - ) (x) in den- 
sity profiles, which varies over the classical Debye screening length, can be 
seen as being enforced by the interplay between the global surfacic elec- 
troneutrality condition and the fact that, when wave-functions vanish over 
different length scales, a charge-density profile appears in the very vicinity 
of the wall even in the zero-coupling limit (see H2.29JI 1, 

The previous calculation can also be interpreted as follows. We no- 
tice that (j qu( - KA ^ and cr qu ^ A ^ can be viewed as the leading fr-terms in the 



38 



contributions to c qu from the regions x < I and x > I respectively, with 
A < I < £ D . Indeed, 



cr 



qu(«:A) 



lim lim / dx > e a \p a fx) 
(A/0-* (VSc)-i-o W ^ L 



(«A) 



(,:) 



(6.15) 



whereas 



qu(ftA) 



lim 



(A/O-o (z/€ 



lim / > e a \p a (x) 



(«A) 



cl(e) 



(x) 



(6.16) 

where p a {x) — p c ^ E \x) is given in H6.12J1 . Z can be identified with the mean 
interparticle distance a, by virtue of (|1.8I) . Therefore, since A <C a, cr q ; u(KA ) 
can be seen as a surface charge located in the plane x = 0, whereas cr 5 u ^ A - 1 
is spread in the fluid over the scale £ D . On the other hand, according to 
1(09)1 and lf6~T3f . $q u ( KA )(a;) may be written as 



qu(«A) ( ' \ _ 



(x) = / dr'a^ (KA) <5(x')0 (O) (r,r'). 



(6.17) 



The interpretation of the latter equation is that $ qu ( KA ) (x) is the classically- 
screened electrostatic potential created by the part of the fluid charge- 
density profile which is concentrated near the wall. 

In the case of an intrinsic semiconductor near a junction, the system 
of electrons and positive holes in the conduction band can be considered 
as a two-component Coulomb fluid of charges — q e and +q e embedded in a 
medium of relative dielectric constant e m . q e is the absolute value of the 
electron charge and energy terms involve e = q e / \ft~^ (see the comment 
after Ijl.llh . Since the semiconductor is intrinsic, the densities p- and p+ 
are equal to each other. They are determined from the energy gap Eq and 
from the effective masses m e 3 and by ^3] 



P = P±(fi) = - A 



3/2 



(6.18) 



Since the system is charge-symmetric, there is no classical contribution to 
the potential drop $(x) and ll.lOt becomes 



p±(x) = p 1 



,-2x 2 /A 



■■■ /\2_ 



<le 



1 - eL{k d x) t /3^$ qu (0)e- KD 



where e = I3q 2 /{2e va ) and, according to 115, 43|) and H5.44I) . 

l /¥ 



1 



(6.19) 



(6.20) 
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(We notice that on principle the expression Ij6.19ll is valid only when the 
wall has the same dielectric constant as the medium where charge carriers 
move.) 

In the case of GaSb, E G = 0.67 eV at 273 if, ra c f = 0.047 m e , mf = 
0.5 m e (where m e is the electron mass), and e m = 15. This system is in 
the regime ill,9Jl for which explicit analytical expressions are calculated in 
the present paper: k d A_ = 2.5 10~ 3 , e = 6 10" 4 , (A_/a) 3 = 2 10~ 6 (with 
(4/3)7ra_p_ = 1). (The length scales are A+ ~ 2.5 nm, A_ ~ 8.3 nm, 
a ~ 640 nm and £ D ~ 3300 nm.) The vanishing of particle densities occurs 
on two different scales A+ and A_ ~ 3.3A + . According to I|fi.l3|l . where 
e a = ±q e , the resultant surface charge located on the wall (over the width 
A_) is cr^ u(KA) = 5 10- 14 Ccm- 2 = 3 10 5 g e cm" 2 . The bulk density of charge 
carriers is p_ + p+ ~ 2 10 12 cm~ 3 and the charge density at distances x > a 
is p c exp[— k d x] with, according to l|6.19Jl and l|6.20ll . p c — 10 9 q e cm -3 . 
Thus p c A- ~ 1 q e cm -2 is indeed negligible compared with o-^. u ' kA ' . The 
potential drop <£> qu (0) = 1.310 _5 ey remains negligible with respect to the 
energy gap. 



7 Comment 

In this section, we comment on the case where the wall is made of a dielec- 
tric material, characterized by a relative dielectric constant e w with respect 
to the vacuum, when e w is different from the relative dielectric constant e m 
of the medium where charges move. Then the Coulomb interaction reads 

v w {r,r') = - A e i, l — (7.1) 

|r — r'l |r — r'*| 

with A c i = (e w — e m )/(e w + e m ). The response of the wall induced by 
the presence of a particle with charge e a (which includes a factor 1 / ^/e^ 
in interaction terms) is equivalent to the presence of an image charge at 
position r*, symmetric of the real-particle position r with respect to the 
wall, and which carries a charge — A c ie a . The Hamiltonian also involves 
the self-energy — A e i e 2 /Ax, due to the interaction of a particle with its own 
image charge. The corresponding loop self-energy can be incorporated in 
the loop fugacity z{C), which now reads 



z(C) = z a d(x) exp 



r ds . ' 



4 J x + \ a £ x (s) 



(7.2) 



In order to exhibit the screening of the self-energy, which is not inte- 
grable at large distances x from the wall, we have performed a resummation 
in two steps, which generalizes the method devised for classical systems in 



40 



Ref.[5J. The choice of the same nine resummed bonds as those in Section 
13.21 combined with the two-step resummation leads to resummed weights 
which are integrable, because they involve only screened loop-fugacities. 

Indeed, in the one-step resummation of Section 13.21 the weights corre- 
sponding to the nine resummed bonds are z sc (C) and z sc (C) — z(£), instead 
of the weights z(C) and z sc (C) — z(C) that arise when there are only five 
resummed bonds without any "double" bond, such as [F cc ] 2 /2, as it is 
done in Ref.[3]. At the end of the two-step resummation process, the con- 
struction rules for resummed diagrams are the same ones as in the one-step 
resummation of Section 3.2, with the only difference that weights z sc (£) 
and z sc (C) — z(C) are replaced by weights z sc ^- 2 \C) and z sc ^ (£) — z sc W (£). 
These weights are integrable at large distances x from the wall, because 
both z sc ^ (£) and z sc M (£) result from resummations of Coulomb ring sub- 
diagrams. 

More precisely, the expressions of the resummed loop fugacities z sc W (£) 
(with i — 1,2) are given by l|3.11|l and i!3,12|l where the value of z(C) 
is that given in Q7.2J1 and <p — v in V^ oud (C) is replaced by cftw — v w ; 
the (pw's have the same boundary conditions as v w written in 117. If) and 
they obey the inhomogeneous Debye equation lll,lfi|) where kM (x) = 
47r/?£ a e£J"2> s ,a(€M£) and ^ (x) = 4n/3 £ Q e 2 a J T> x , a (t) z sc W (£). 
(Resummed bonds are defined with (p l w ' in place of (j>.) 

As a consequence of quantum dynamics, screened loop self-energies (and 
subsequently particle densities) are found to approach their bulk values 
only with an integrable l/x 3 tail, whereas the particle self-energy due to 
the electrostatic response of the wall is exponentially screened in classical 
systems ^HQ]. Even for bulk properties 53Ej, screening in quantum 
systems is less efficient than in classical fluids. 

More precisely, the screened loop self-energy in z sc W (£) or z sc ^ (£) is 
the sum of two contributions. The exponentially-decaying part has the 
same decay at large distances from the wall as the screened self-energy of 
a classical charge 0. The algebraic l/x 3 tail arises from the other part, 
which reads 

A 4r^4-*<'--'0 _w^__ ( , 3) 

(See the analogous term l|3.9|l in the screened pair interaction). In the low- 
degeneracy and weak-coupling regime the particle density p a {x) does not 
seem to have a simple explicit value, because of the self-energy contributions 
arising from the dielectric response of the wall. 
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A : Resummation of large-distance Coulomb 
divergences 

In this appendix, we display the resummation of Coulomb divergences in 
the Mayer representation (|3.2I) for the loop-fugacity expansion of the loop 
density. We use the same decomposition of /-bonds as the one performed 
in Appendix B of Ref. 0] . (The a priori arbitrary decomposition is chosen 
according to the properties to be studied after resummation.) The /-bond 
is written as the sum of ten auxiliary bonds /, 

/(£,£') = 

+\ [f cc f + / c 7 mc + / c 7 cm + \ [/ mc ] 2 + \ [/ cm ] 2 + / TT | (A £')• 

(A.i) 

Bonds / cc , / mc , / cm and / mm are defined from the multipolar decom- 
position {HUH) and are introduced in order to handle classical exponential 
screening. These bonds read 

r b (£, £') = -f3e a e a , V ab (£, £'), (A.2) 

where superscripts a and b stand either for c or m, and where we have 
set V cc (£,£') = v(r — r'),. / TT is defined by (|A . 11 . Double /-bonds 
(such as [/ cc ] 2 /2, / cc / cm , ...) are not involved in the series representation 
!■'•>. 2t where two points can be linked by at most one /-bond. We choose 
to make them appear in the decomposition llA.l(l in order to obtain a 
finite sum for the contribution from all possible ring subdiagrams which 
are defined hereafter. When the decomposition l|A.lj) is introduced in the 
Mayer representation ll.'i.2t of the loop density, diagrams G are replaced 
by diagrams G built with the ten /-bonds defined in <| A . 2ll and the same 
topological rules as diagrams G. 

The resummation procedure relies on the integration over the inter- 
mediate points (called "Coulomb points") of all possible Coulomb-chain or 
Coulomb-ring subdiagrams. A Coulomb-chain subdiagram between two 
points £ and £' is equal to 

J d£x ■ ■ ■ d£ N f ac (£, A) z{£x) f cc (£ u £ 2 ) z{£ 2 ) f cc {£ 2 ,£ 3 ) 

x ■■■z(£ N )f ch (£ N ,£') (A.3) 

with an arbitrary number N > 1 of internal points. A Coulomb-ring sub- 
diagram attached to a point £ is a closed Coulomb-chain subdiagram, i.e. 
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+ f 



with £ = £' (see Fig®. It is equal either to / d£ x z{£ x ) [/ cc (£, A)] 2 A 

/rf£ l2 (A) [/ mc (£,A)] 2 /2 , /^i2(A)/ cc (/:,A)r c (A4) or 



where 5 is the symmetry factor of the Coulomb-ring subdiagram (S = 2 
if a = b and 5=1 otherwise). Coulomb- ring subdiagrams do exist in 
diagrams G, because the latter ones contain articulation points. 

Diagrams G can be collected into classes of a partition where all di- 
agrams G inside the same class lead to the same prototype diagram P* 
after erasing the intermediate points of all Coulomb-chain or Coulomb- 
ring subdiagrams. When resummed bonds are defined, some "excluded- 
composition" rules ensure a one-to-one correspondence between each class 
in the partition of diagrams G and each diagram P*. In the present paper, 
we choose to build prototype diagrams P* made with bonds F cc , F cm , 
F mc , (1/2) [F cc ] 2 , (1/2) [F C1J1 ] 2 , (1/2) [F mc ] 2 , F cc F cm , F mc F cc and F RT 
(see Q3.7p ). This choice leads to resummed diagrams where the new weight 
of every point is convenient for dealing with the case of a wall with a 
dielectric response, as explained in Section 

The bond F cc (£,£') between two loops £ and £' of a prototype dia- 
gram P* is defined as the sum of the bond f cc (£, £') and of all Coulomb 
chains involving only / cc -bonds between £ and £' (see Fig. ®. The bond 
F mc (£,£') is the sum of the bond f mc (£,£') and of all Coulomb chains 
fOjl with a = m and b = c (see Fig®. Bonds [F cc ] 2 /2, [F cm ] 2 /2, 
[F m f/2, pccpcm and pccpmc between tw0 points £ and £ originate 

from subdiagrams with either one bond (1/2) [/ cc ] 2 , (1/2) [/ cm ] 2 , 
(1/2) [/ mc ] 2 , / cc / cm , or f cc f mc respectively, or from subdiagrams with 
two Coulomb chains between £ and £' (see FigQJ. In Ref.[3|, the proto- 
type diagrams were chosen to be built with only five resummed bonds called 
pec ^ pem^ pmc an£ j , Since the resummation procedure relies 

on the integration over the same intermediate points in the same diagrams 
G in both cases, the expressions of F cc and F mc are the same (see Eqs. 
H3.7ajl and Il3.7hjl ). while F RT is equal to 



F RT = F Rz --[F cc ] 2 --[F cm } 2 --[F mc ] 2 -F cc .F cm -F cc .F mc . (A.5) 



>From the expression of F Ri derived in Ref. 0J, we get the value l|3.7cll . 

The "excluded-composition" rule is different from those of Ref.@], be- 
cause resummed bonds are different in the present case. The rule arises 
from the following arguments. When £ is the intermediate point of the 




dCi • • • d£ N / ac (£, A) z(£x) f cc (£i,£2) z{£ 2 ) f cc (£ 2 ,£ 3 ) 



x ■■■z(£ N )f ch (£ N ,£) (A.4) 
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chain F ac (d,£) F ch (£, Cj) with i ^ j and is not involved in any other 
bond, or when £ appears only in one bond (1/2) [F ac (d, £)] 2 (with a = c 
or m) or one bond F cc (d, £)F mc (d, £), the reason why £ has not dis- 
appeared in the resummation process can only be that £ carried at least 
one Coulomb ring in every diagram G in the class which leads to the di- 
agram P* after erasing Coulomb points. Therefore, after integration over 
Coulomb points, the weight z{£) attached to point £ in diagrams G is 
multiplied in diagrams P* by the sum of all products of Coulomb rings 
attached to £. As argued in Ref.jlj, the sum is equal to expj r — 1, where 
I r is the sum of all possible Coulomb rings (see Fig. EJ, The corresponding 
weight is equal to z sc (£) — z(£) where z sc (£) = z(£)exp/ r . I r is linked 
to the screened potential <f> by I T = — Pe^V^ oud where V^ c oud is defined in 
H3.12I) . z sc (£) exhibits the stabilizing effect of screening clouds t 5.j. When 
£ is not involved only in a product F ac (d,£)F ch (£ 7 £j) (where d may 
coincides with £j ) , £ may carry no Coulomb ring in the diagrams G of the 
corresponding class, and its weight in diagrams P* is equal to z sc {£). 



B : Scaling analysis 

In this appendix, we perform the scaling analysis of diagrams with one 
internal point in the resummed Mayer expansion H3.10I) . 

B.l F mc - and F cm -bonds 

A simple scaling analysis of the contribution from diagrams with either one 
bond F cvn (£,£') or one bond F mc (£, £') is obtained by a Taylor expansion 
of their expressions (see l|3.7bjl l around their classical values. 

In the case of diagram with one bond F m c the Taylor expansion in the 
variable A a £(s) can be factorized: the leading order is given by 

J dsA Q |(s).V r (J d£'F cc {£,£') z sc (£')j ■ (B.l) 

According to the discussion in Section |Ol at leading order the parenthesis 
proves to be a function of r which varies over the scale k _1 . By using the 
fact that the action of the gradient A a V r on a function of ky multiplies its 
order by k\ , we get from 115. 2 2|) that 

J d£'F mc (£,£') z sc {£') = o[e.n\{n\) 2 y (B.2) 

The diagram with one bond F cm (£,£') involves integrated moments 
of £ , which depend on A 7 , and e 7 z 7 can no longer be factorized at 
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leading order. Moreover, the first moment of the ^-component of £ does 
not vanish, and the order of J dC F cm (£, £') z sc (£') is given by 

5> 7 | dr' Q\s J 2WO&(*)) A / FC ^' r,) (B.3) 

At leading order X 7 dF cc (r,r')/dx' is a function of nr and fix', while, 
as can be checked in i!2,18|) . the mean extension of a Brownian path, 
Jq ds J D x /^(£')£' x ,(s), vanishes gaussianly fast when x'/Xj goes to in- 
finity. The latter function makes the next x'-integration convergent over 
the scale A (and not k _1 ). Therefore, as in l|5.23[l . the order of the integral 
ID .31) is equal to kA times the order of z 1 J dr' A 7 d.F cc (r, r')/dx'. Since at 
leading order F cc (r, r') is a function of nr and nr', the latter integral is of 
order kA times the order of z 7 J dr' F CC (C, £'), which is equal to 0(1) by 
virtue of 115. 21|) . Eventually 

J d£'F cni (£, £') z sc (£') = O ((kA) 2 ) . (B.4) 



B.2 Multiple bonds 

In the case of multiple bonds, the dependence upon the species of the loop 
£ is not reduced to e 7 z 7 . Therefore, contrary to what happens in the case 
of F cc , the neutrality constraint on fugacities 12.1111 no longer increases 
the actual order with respect to that given by pure scaling analysis. For 

instance, since [F cc ] is proportional to e 2 </>(°) (?,?') 

portional to e</>(°)(r, r'), H5.21H implies that 



, while F cc is pro- 



J dC i [F cc ] 2 (£, C')z sc (C) = 0(e). (B.5) 

However, owing to the construction rule about fugacities, the internal 

point in the diagram made with one bond [F CC Y /2 has a weight equal to 
[z sc (C) — z(C')]. By virtue of H5,15[l . the latter weight is of order O(ze) 
and l|B.5p leads to 

J dC! 1 - [F cc ] 2 (£, C) [z sc (C) - *(£')] = 0(e 2 )- (B.6) 

Similarly diagrams [F cc F mc ] (£,£') and \ [F mc ] 2 {£,£') also have an 
internal point which carries the weight z sc {£')~ z(£'). Only the moments of 
the Brownian path £, which is associated with the nonintegrated position 
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x, are involved and, according to the argument leading to llB.2li and by 
virtue of lIB.fij) . the contributions from these diagrams are 

and 

J dC'^ [F mc ] 2 (£, £') [z sc {C) - z{£)] = O (e 2 . (kA) 2 ) . (B.8) 

For the bond [F cc F cm ], by virtue of the mechanism involved in l|B.4ll . 
the integration over x' of the first moment of the Brownian path £' leads 
to an extra factor (kA) 2 with respect to the order JB.51 of the contribution 
from [F cc ] /2, namely 

J dC [F cc F cm ] (£, £') z sc (C) = o(e. (kA) 2 ) . (B.9) 

At leading order the bond [F cm ] 2 /2 involves the second moment of the 
Brownian path This moment tends to a non-zero value when x goes to 
infinity, so that the integration over x' still converges over the scale 
Therefore, according to the argument leading to JB.2L the order of the 
contribution from [F cm ] /2 has an extra factor of order (kA) 2 with respect 
to ilB~5j) . 

J dC'^ [F cn f (£, C) z sc (C) = 0(e. (kA) 2 ) . (B.10) 
B.3 F RT -bond 

Let us consider the diagram made of a single bond F RT . The integration 
over the Brownian path £' makes the integral over r' convergent at small 
distances |r — r'|. The order of the contribution from this diagram can be 
inferred from the known results about the bond F* T obtained by adding 
to F RT the double bonds different from [F cc ] 2 /2, which do not contribute 
to the integrability of F RT at large distances |r — r'| in the limit where e 
vanishes. F* T is defined by 

F RT = F* T - - [F cm ] 2 - - [F mc ] 2 - F cc .F cm - F cc .F mc . (B.ll) 

where F RT is given in l|3.7c|) . 

The order of the contribution of / dr' J T> XiCt (£) J O x ', 7 (0 F* T (£,£') 
has been studied in the bulk situation in Ref.|19j (and a similar calculation 
also appears in Ref.[20j). Let us introduce the thermal de Broglie wave- 
length X ai associated with the reduced mass m a m 7 /(m a + m 7 ). In the 
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bulk, the considered integral proves to be the sum of a term of order A 3 
times an analytic function Q(£ Q7 ) of the parameter £ a7 = — /3e Q e 7 /A Q7 
[Tl] . plus a term of order (/3e Q e 7 )A 2 (with two contributions where A = X a 
or A = A 7 ). The latter "diffraction" term arises from the second moments 
of the Brownian paths. There is no term proportional to (/3e Q e 7 ) 2 A, be- 
cause the first moment of a Brownian path vanishes in the bulk, by virtue 
of the spherical symmetry. 

In the vicinity of the wall, the precise calculation is more delicate than 
in the bulk. However, we may expect to obtain the same orders AjL, and 
(/3e a e 7 )A 2 as in the bulk, plus a term of order (/3e a e 7 ) 2 A allowed by the 
anisotropy introduced by the presence of the wall. Since the fugacity z sc (C) 
is of order O(p) = 0(a~ 3 ), the contribution from the diagram made of a 
single F RT -bond can be viewed as the sum of four terms with respective 
orders 

A \ ^ i /3ea£-t 



Gi / \ '"CK 7 



~ Qw , (B.12) 



where Q w is defined similarly to Q mentioned previously pi] with the only 
difference that space integrals are restricted to x > and x' > 0, 



0(( K A) 2 ), (B.13) 



a I \ a 



- ) =0(e.k\), (B.14) 



and 

' R x 3 

fj 6q G/y 



ln(^A) = O (e 2 Iii(kA)) . (B.15) 



a 

In these equalities we have used the relations jl.fijl and |QJ . 



References 

[1] J.-N. Aqua and F. Cornu. Density profile in a classical Coulomb fluid 
near a dielectric wall. I. J. Stat. Phys., 105:211, 2001. 

[2] B. Jancovici. Sum rules for inhomogeneous Coulomb fluids, and ideal 
conductor boundary conditions. J. Physique, 47:389-392, 1986. 

[3] A. Alastuey, F. Cornu, and A. Perez. Virial expansions for quantum 
plasmas: Diagrammatic resummations. Phys. Rev. E, 49:1077, 1994. 

[4] F. Cornu. Correlations in quantum plasmas. I. Resummations in 
Mayer-like diagrammatics. Phys. Rev. E, 53:4562, 1996. 



47 



[5] J.-N. Aqua and F. Cornu. Density profile in a classical Coulomb fluid 
near a dielectric wall. II. J. Stat. Phys., 105:245, 2001. 

[6] E. Lieb and J.L. Lebowitz. The constitution of matter : Existence of 
thermodynamics for systems composed of electrons and nuclei. Adv. 
Math., 9:316, 1972. 

[7] F. Cornu. In preparation. 

[8] M. Kac. Probability and Related Topics in Physical Science. Wiley 
Interscience, New York, 1959. 

[9] B. Simon. Functionnal Integration and Quantum Physics. Academic, 
New York, 1969. 

[10] J.-N. Aqua and F. Cornu. Classical and quantum algebraic screening 
in a Coulomb plasma near a wall : a solvable model. J. Stat. Phys., 
97, 1999. 

[11] E. Meeron. Theory of potentials of average force and radial distribu- 
tion functions in ionic solutions. J. Chem. Phys., 28:630, 1958. 

[12] J.-N. Aqua and F. Cornu. Dipolar effective interaction in a fluid of 
charged spheres near a dielectric plate. Phys. Rev. E, 68:026133, 2003. 

[13] D. Zwillinger. Handbook of differential equations. Academic Press, 
1989. 

[14] W. Ebeling. Coulomb interaction and ionization equilibrium in par- 
tially ionized plasmas. Physica, 43:293, 1969. 

[15] Ashcroft N.W. and N.D. Mermin. Solid State Physics. Holt-Saunders 
International, 1976. 

[16] L. Onsager and N.T. Samaras. The surface tension of Debye-Hiickel 
electrolytes. J. Chem. Phys., 2:528-536, 1934. 

[17] A. Alastuey and Ph. A. Martin. Absence of exponential clustering in 
quantum Coulomb fluids. Phys. Rev. A, 40:6485, 1989. 

[18] F. Cornu. Correlations in quantum plasma. II. Algebraic tails. Phys. 
Rev. E, 53:4595, 1996. 

[19] F. Cornu. Quantum plasmas with or without a uniform field. II. Exact 
low-density free energy. Phys. Rev. E, 58:5293, 1998. 

[20] A. Alastuey, F. Cornu, and A. Perez. Virial expansions for quantum 
plasmas: Maxwell-Boltzmann statistics. Phys. Rev. E, 51:1725, 1995. 



48 



